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ABSTRACT 



The theory of the three-dimensional rotational flow of an in- 
compressible and inviscid fluid through an axial turboraachine is 
described end the hydro dynamical equations are simplified by consid- 
ering an infinite number of blades in each row. The forces of the 
blades on the fluid are treated as non-conservative body forces dis- 
tributed uniformly about the axis. 

Formulation of the mathematical problem leads to one non-linear 
partial differential equation and two integral equations for the 
throe velocity components. A linearized solution of those simultan- 
eous equations for any prescribed blade loading is based on the 
consideration that the vorticity gener ted by the blades is trans- 
orted downstream by the mean axial velocity. An iteration process 
which leads to solutions of greater accuracy is developed by consid- 
ering for each iteration that the vorticity is transported by the 
velocities found "by the previous iteration. 

The Bessel’s functions w..ich occur in the Green’s function solu- 
tion are replaced by their asymptotic values and the infinite series 
is summed to express the solution in closed fora. The iteration 
process is then adapted to mechanical calculations by dividing the 
region of vorticity into small rings of rectangular cross-section and 
determining the influence on the velocity of a unit change of vortic- 

i 

ity in e ch of these rings. Once this influence is established it 
is relatively easy to calculate the velocities in any axial flow 
machine with any prescribed blade loading. 
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INTRODUCTION 



The development of turbomachinery until recent years has resultod 
primarily from empirical methods, both in design and analysis. Recently 
theoretical methods have led to the advancement of axial flow compressors 
and turbines even though, because of the complicated mathematical problem, 
□any simplifying assumptions were necessarily made and solutions were 
approximate. 

In analogy with the three dimensional wing theory, problems of 
flow in turbomachines may be classified as: 

(1) The Direct Problem ! The direct problem of calculating the flow 
In turbomachinery is that of detexmining the velocity field, 
the blade forces, and the distribution of energy in the fluid 
when the blade shape, the blade speed, end the appropriate 
boundary conditions are prescribed, 

( 2 ) The Inverse Troblenu The inverse problem of calculating the 
flow in turbomachinery is that of determining the velocity 
field, the blade shape, and the distribution of energy in the 
fluid when the blade loading, the blade speed, and the boundary 
conditions are prescribed, 

(3) The Intermediate Problem An intermediate problem which ap- 
pears to be of interest is that of determining the velocity 
field, the magnitude of the remaining forces, and the distrib- 
ution of energy in the fluid when one blade force is prescribed 
and the blade shape is partially prescribed. 



2 - 



The diroot problem arises when it is desired to investigate a 
given machine "off the design point*. 

The inverse problem nay bo solved for the Initial design but 
usually, because of the structural limitation on blade shape (see 
appendix), it will be more practical to formulate and solve the inter- 
mediate problem. 

In order to simplify the mathematical problem the fluid will be 
assumed to be inviscid and incompressible and the blade forces will 
bo treated as body forces, uniformly distributed through the fluid, 
so that the flow is symmetrical about the axis of rotation and the 
vorticity is no longer shed in sheets behind each blade but is contin- 

ft 

uou3ly distributed over tho region downstream of the blade row. 

An axial flow machine for which the inner and outer boundaries con- 
sist of concentric circular cylinders extending to infinity in the di- 
rection of the flow will bo considered. 

Tho mathematical problem is formulated by considering the time 
rate of change, along a stroamsurfece, of tho tangential component of 
the vorticity vector. The diffic dty of this problem lies in the 
solution of the non-linear partial differential equations that describe 
rotational fluid motion. In order to overcome this difficulty a method 
of iteration is developed whereby solutions of any required degree of 
accuracy may be obtained. 

The first step of this iteration process provides a linearized 
solution (cf. Marble*) based on the assumption, analogous to the 
Prandtl three dimensional wing theory, that the vortioity is "trans- 
ported* downstream by the mean axial velocity and is not influenced 



by its own induosd velocities. This linearized solution is not as ac- 
curate as might be expected from the above analogy, one reason being 
that the vorticity is shed in three dimensional space instead of in 
a two dimensional sheet so that the induced velocities are likely to 
be veiy large. 

for the seoond step of the iteration process the vorticity i3 
considered transported by the velocities found by the linearized so- 
lution of the first step. The succeeding approximations are obtained 
in the same way, in each case using the velocities of the preceding 
approximation. 

The solution is obtained by finding the appropriate Green's func- 
tion - that is, a function G (r, z;«<,^ ) which gives the velocity, 
consistent with the boundary conditions. Induced at any point of a 
circle r, z by a unit change In tangential vorticity of a vortex ring 
at radius at an axial coordinate / 3 . The Bessel functions which 
arise in this solution are replaced by their asymptotic values and the 
Green's function, which would occur as an infinite sum of Bessel func- 
tions, is bythis moans expressed in closed form. 

The resulting expression for the radial velocity, a double inte- 
gral, is modified and adopted for mech&nioal calculation. 
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I. Notation and : 
The flow is described (Fig. 1) in 
r, , z, by the velocity components 
u, v, w, respectively. The corres- 
ponding radial, tangential, and axial 
velocity components are 

n v 

S “ o z. 



lymbolc 

a cylindrical coordinate system 

u,S 





D w 

O r 



5 



j_ 1L_ fro) 
^ 5r 



Fig. 1 

Coordinate Gy stem and Desig- 
nation of 'Velocity and Vortio- 

ity Components 



In addition ?;e vili uso the following symbols. 



cO - angular velocity of rotor 
p - pressure 

V r absolute volocity vector 
= vortioity vector 

F 1 force vector (foroe ox' blades on fluid) 
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II. Ilydrodynamic Equations 

The following equations are simplified for this case of an invisoid 
and incompressible fluid in steady, adiabatic, axially symmetrical flow. 



Equations of Motion : 



V ■ v V 
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7 - if - f 



( 2 ) 



Continuity Equation i 



p • 1/ - o 



(3) 
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III. Formulation of the Mathematical Problem 



a) Development of the Basic Zfauatlono 

Because of the axial symmetry only the tangential vorticity is 
associated with the radial end axial velooities ( while the radial and 
axial vorticity components ere only associated with the tangential 
velocity* The tangential vorticity constitutes an annular vortex ring 
and by considering the deformation of this ring information concerning 
the radial end axial velocities oan be gained* It is well known (of* 

p 1 

Meyer , Marble ) that the circulation about a deforming vortex ring ie 
oonotant in a conservative force field. Only a radial stretching oan 
occur sinoe the flow is oymaotrlcal about the axis so that the constancy 
of circulation requires that in a conservative field 



In the presence of non-oonservative blade forces or if radial 
vorticity ia present the circulation, and henoe the quantity will 
not remain constant. The law governing its variation oan be derived 
from the equations of motion. The equations of motion in veotor form 
can be written 



If ve take the curl of both sides and use the equation of continuity. 




( 4 ) 



v * 12 = yr * f ■ r - p 



( 2 ) 



cr . 1 / 



o 



( 3 ) 
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there re suit a 

/ • ? Q - r \/ ~+“ v + F 



The tangential conponont of this vector equation is 



u. 



Djl 

Dr 



2n 
d z 




Dv 



r ^ 



CL 



df<r _ D^z 



which, using the definitions of vorticity, can he simplified to 

£,(*)- i-fn (? )*&-&/ . , < 5 > 

Since the circulation i3 directly proportional to the quantity-- , 
this equation expresses the law governing the tine change of circula- 
tion around the tangential vorticity ring. The first term on the right 
io the change in the axial direction of the centrifugal force, the sec- 
ond term is the change of the radial blade force in the axial direction, 
and the third tern is the change of the axial blade force in the radial 
direction. A little thought will show that these terns, in each case, 
represent monents tending to cause rot-tion of a particle about a tan- 
gential axis. It is in thi3 manner that a non-conservative force field 
tonds to effect a change in the tangential vortex ring. 

The tine derivative on the left is taken along a streamline and 
is written 



" k(Z) 



( 6 ) 
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In order to compare the relative magnitude of the terms on the right 
we use the definition of tangential velooity and write 






&<2) = 



w 



d z 



(II 



K t) « _ U D *4/ _ C*- Jj£ y. U. 7) vv 

/*”* 3772. r~ 3 /-£ /~z- 3 i. f-z & r~ 



M/ c) x a _ \ a / D 

r Dz*- r Ttjz 



Remembering that for the axial flow machine the radial velooity is 
small compared to the axial velocity and expecting from physical con- 
siderations that the velooity distribution will be smooth, it appears 
that 

u £- r (-£} < ^ " il(-r) 



This inequality will be useful for the first approximations but must be 
more closely investigated for the final approximations. It is conceiv- 
able that for flow that differs greatly from vortex flow the radial 
change of tan<~ontial rorticity would be of such magnitude that this 
inequality is not Justified. It will however be useful to group the 
smaller terrrn separately, that is, to write 



DZ 




Dr 
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using the definition of vortioity we can write an equation for the 
radial velooity as follows 



27 * 




D^u — _L 
+ 31* ~ ^ 
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Dft 

Dz 



J Fz _ u 2 . _l( dtt. 7 

Or dr r { az a r / J 
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ao that the left side, if set equal to zero, would be a linear partial 
differential equation for which a solution is known. 

Since the fluid is inviscid the force of the blade on the fluid 
will act normal to the blade and hence normal to the relative velocity 
so that 



U /> F+ 



T 4 * PV' 



F? = ° 



( 9 ) 



The radial velooity and the radial force are both small compared to the 
other force and velooity components so that with very good approxima- 
tion 



1/ - ^ r — 

Fz = F-& 



\AS 



( 10 ) 



Using thl s relation in the equations for radial velocity we have one of 
the final equations of the iteration process 



cFu. 

Pr*- 



/• 



D fU I Du 
sA'rJ+V#. 




d Fr D / / Du _ 1 1 

3Z rl'fl d r )J 



(ii) 



For determining a relation between the tangential velooity and the 
tangential force we havo the second oquation of notion 

(rr) = +"§z( r '') = r F^. 

cf r 

which states that the variation along a streamline of the moment of 
momentum is equal to the mcnent of the tangential farce about the axis 
of rotation. Then since we will prescribe the tangential force we 
can express the tangential velocity as an integral along a oylinder. 
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u. 





( 12 ) 



Th© continuity equation. when integrated# provides an equation for 
the axial velocity in terns of the radial velocity 






I AJ 
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'a - - 



/ 
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2-(ru) 

dr 



c//3 



( 13 ) 



V7e have developed three basic equations, Eqa. 11, 12, 13, to be 
used in formulating the linearized problem, for the first approximation, 
and in the eonatructlon of an iteration procedure whereby more exaot 
approximations nay bo obtained. 
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b) the Linearized Problem 

In order to linearize the basic equations derived in the previous 
section we shall consider axial flow, where the radial velooity is very 
small compared to the axial velocity, 

u < < w 

and will assume that the vortioity is transmitted axially downstream 
with the mean axial velocity w o so that 

(~r ) ~ "K* §z ( 2) 

Furthermore, if the radial force is small, we have for calculating the 

radial velooity the expression 

2 *~^t- 1 D 1 u 

3r * <£> r ‘ J 3 z. 2 - 

+ U4) 



For calculating the linearized tangential velocity Eq. 12 becomes 



^ z f s>- ( r- y 

v = + J ~1 aT 0 ^ P 

/S = -•*» 



and the linearized axial velocity is, from Eq. 13* 



U5) 



~ - / yr h (ru) ^ ( 16 ) 

73 = - 00 

We are now able to outline a solution to the linoarized problem. 

We will consider here the inverse problem and will specify the two 
force components F r and . *'e will prescribe F r <.<■ F^. and F^- F^Xr.z), 

«r>d will specify the boundary conditions applicable to axial flow. 
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First the tangential velocity is determined by 



V 



t'.. 



r 

/ 
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showing the tangential velocity we can determine the radial velooity 



from I Iq. 14 



d 2 j± 

s r *• 



- /V,z; 

— / A 2~ (jl~ 1 , 

kV, l_ 2 z ' r *♦'# 
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with the boundary condition that 



u = O 



/■ = 






r t 



z = ^ 



We have here an olllptlc partial differential equation for the radial 
velocity. Since we know the right hand aide, the unhomogeneous part, 
and have complete boundary values we can solve this equation for the 
radial velocity. This equation was solved by Marble^ in terms of a 
Green's function G(r,z{ oc,/? ) which gives the velooity induoed at any 
point of a circle r,z by a unit chance in tangential vortloity of a 
vortex ring of radius ^ at an axial coordinate /? . By Marble's solu- 
tion the radial velocity is expressed as 



a. 




G (r, z • ^ /3) c/~c a/ 0 



(16) 



where 



/?<*, C3) = ^ [r?( r ) + 7777 r ^a-] j ^ 



(17) 



G(c*;«p) = ^ 



n - / 



a fa*) e 

2 



-«-„/z-W 



Z 



( 18 ) 
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U^(e n r) is the linear combination of Beasel functions of order one 



LT,{£ n r) ~ J, (±*r) )^ { e n r,) - J, (€■*/■ J y,(er„r) 



( 19 ) 



and the characteristic values e„ are chosen for those particular bound- 
ary conditions. The quantity ^ is the nom of U^fe/) over the in- 
terval rj_, r 2 . 



\ 2 2 ( 6 n /a ) ~ K (^n 

= £ 



(30) 



Knowing now the radial and tangential velocities we have from 

Bq. 16 the linearised axial velocity 

(■ 3 - z 

-jr jp ('“I dp 06 ) 

r _ c^» 

We have formulated above the linearised solution as developed by 
Marble 1 . It should be noted that the Green's function as derived is 
independent of the manner of linearization provided, of course, that 
the unhonogeneous part f(<*,0) is known. If then we have another value 
of f (~<, P ), more exact than the above linearized value, we can use the 
sane integral (Bq. 16) to determine a more exact radial velocity dis- 
tribution. 

This integral, involving an infinite sum of 3essel functions, is 
extremely difficult to calculate. A solution for a particular blade 
loading was obtained by Marble 1 . For the axial flow machine where the 
boundaries exclude the region near the axis, l.e., the small values of 
the arguments of the Bessel function, we con use the asymptotic values 
of the Bessel function with good accuracy. This makes it possible to 
sum the infinite series so that the Green's function can be expressed 



in closed form. 
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c) Asymptotic Value of Groo n 'a Function 

An asymptotic approximation, in closed formula desired for the 
Green's function 



£ ( *■ i 

where 




-e n /z-t 9* 

V, (e„r) 

2. ± n *V 



( 13 ) 



U,((r n r) - J, (6 n r) / t (e„ rj - J, (e n r.) K (e„r) 
V a (£„>-) - J Q (£»r) r o fa n r,) -J 0 (e»r.)y 0 (6„r) 



( 19 ) 



T ^ _ cL (to. r *)- (±1 r A (20) 

£ 

and Id found from tho boundary oondition that 



V, (e n r 2 ) = J r (e n r z ) Y,(&„r,) - J, fe n r t ) K C*»rJ = O (21) 



For "large* (e„r) asymptotic valuoa of tho Dosool functions are 
found in Jahnke-3nde^. 



c/, r) 




, 377 . 

cos (€„r - + ) 

Jg- rre n r 




c/o f±n r] 




cos (e n r - 

j£ rre„r 


( 22 ) 


r, C± n r) 




Sin (€ n r- jr) 

J± rre n r 




K (*n H 




s/n (&» /* - 

/Yrre„r 
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Tho co responding characteristic values t„are found from the 
Golution of Sq. 21, with the Bessel functions replaced by their asymp- 
totic values, ::<\ . 22. 

77* 



€/7 n - r , > 






(23) 



Using those oharactoristio values in Uq. 22 wo find the complete asymp- 
totic value of the Bessel functions, compatible with the boundary con- 
ditions of the problem. 



lJ/ (& n f) 

d 0 ( £ n r) 
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//7 rrr 3jr\ 

inlru-r, +■ j 



(24) 
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The asymptotic value of T Jy , and follow os3ily and the asymptotic 

value of the Breen's function oan bo written r.3 the infinite sum 

./z-J*/ 



G = 



'W £ S/ " **'■ 



r- r, 



S/ft ft ft 



r> - / 
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r* -r, 



— r> 7T 

3 



ri-r r 



This series can bo suimed and er.prosocd in closed form. 
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d) Tho Iteration i Procedure 

Using the asynptotio value of tho Green's function a o derived in 
tho previous aoction throe general equations can be written for the 
iteration procedure. 

io will use tho volocitioo resulting from a solution of the lin- 
earised problem to obtain the seoond approximation and will obtain 
subsequent approximations using velocities of the preceding approxima- 
tion. Here we are. in effeot, assuming that the vorticity is, in eaoh 
oeso, transported by the velocities of the preceding approximation. 

The radial and tangential foroe components are prescribed and we 
start with 



»0= o 

v a -- mean axial volooity 

If the velocities which are to be obtalnod by the /7th approximation are 
denoted by tho subscript n the oquationa for tho /rth approximation vel- 
ocities are 
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f3 - Z 




h- dfi 



(28) 



These three equations, 26, 27, 28, are the basis of the iteration 
process which is adapted to mechanical calculations in the next section, 
w'e will make use of the important fact that G(r,s; <*/3 ) is the sane 



for each iteration 
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IV. /Captation of the Iteration Process to 
Mechanical Calculation Methods 

The Integrands of the expressions for the tangential and axial 
velocity components, Dqs. 26 and 28, are regular over the regions con- 
sidered and can be integrated by any numerical method with relatively 
little difficulty. Mechanical calculations can be performed using the 
equations in their present form. 

The integration of the equation for the radial velocity Eq. 27 
becomes rather involved by any method, and for mechanical calcula- 
tions the form of the equation must be radically changed. An inves- 
tigation of the Green's function, 0(r,zjc*,/3) will reveal a logar- 
ithmic singularity at c<= r, p- z. To determine the radial velocity 
at any point r,z the integration must be oarried out over the entire 
region. The logarithmic singularity of the integration represents 
the influence at r,z of an axial change of vortioity at r,z. 

It will be advantageous to perform the integration in three 
parts, one part being over the regular region, away from the singular- 
ity, and the other two parts being over a small region which includes 
the singularity. Integration over the latter region containing the 
singular point consists of evaluating a regular part and evaluating 
a part containing a purely logarithmic singularity, tfith this in 
mind wo write the radial velooity as the sum of three integrals. 

^ u' + u" -t- u"' 



u 



(29) 
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1 ^,—'^ -1 ( 31 , 
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2 ye ry/ 



2-e r-<T 1 ' 



( 32 ) 



Obviously the integrand for u' is finite but wo need to show that 
the integrand of u'' is regular. It is easily shown that 
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7L-1 - 

•'x.-r, 



Jr- rAy 

r, z - r. 



costysr'LJ- 

/l -c /\ - r, 



r&y-m-f 



/•-r. 



— J_ / 77 / / - c-os 2 7T - 7 -, 

4 77 I 



( 33 ) 
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There is a singularity in this integrand, however, but it does not 
occur near = r,z. The distance from r,z to the singularity, oeas- 
urod radially, is given by 



From those equations this singular! ty is seen to lie outside of any 
region whose center is at r,z. The Integration for u" therefore does 
not involve a singularity. 

The integrand for u'" obviously possesses a purely logarithmic 
singularity at /iir,z. 







a) evaluation of the Integral for u' 
Osing notation that 






we have 
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(35) 



The function f («*, 0) will be different for each iteration but the 
Green's function, G(r, 2 j ) will be the same* 
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Fig. 2 



Subdivision of Flow Field 
Into Steal 1 Rectangles 
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If we divide the region into snail reotangles as indicated in 
Fig. 2 we can assume with good acouraoy that f ( <*■, P ) can be repre- 
sented in one of these sufficiently snail regions n by the second 
degree surface 

tC, SV = C„, -a C nz °< C n } ft y 06} 

where four values of f ( **, /? ) in the region n are necessary to eval- 
uate tho constants C. Suppose the four values of ) are 

f . — f , . Substituting those values in Sq. 3& we have four eque- 

n.L TV ^ 

tions from which tho constants C can be evaluated. 
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— C . must be determined for each rectangle. 
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It will bo convenient to use f( ^P) evaluated at the corners of the 
rectangle for this. Using the above notation we can write 



u 



n 




Cm / C , 7 ? *=K + C/?3 /3 + C t 



*7 + 



/&) G (rz- •<! 0) o/a 0//1 ( 3 ^) 



and can further define 



In, = 
JnZ - 

-A 3 ' 

1/7 + 



ff* «„ o'- ^ 

*/-< 

//, /? G„ o'- dp 



(40 
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30 that 



U' = 



{(-/}/ 7m c nZ 7/iZ. + t-nj 7^3 ■/■ C 0 + In* j 
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Ifare the °nl depend only on the function f ( P ) and mu3t be deter- 
mined for each iteration. The I ni depend only on the Green’s function 
and can be determined once and for all. It should be noted however 
that the 1^ are different for each r,z> that is, we must have a cer- 
tain set of I . to use in determining the velocity at each point, 
ni 

m 

If we denote by u* the velocity at the center of themth rectangle 
and use the superscript m to denote the corresponding 1^ % then the 
velocity u* at the center of the /77th rectangle is 



u ' 



n 
n / m 



F _ ' r7 r m j m 7 m 1 

j c n/ J + <-nz In 2 + Cn3 j h * + C'n+J n + j 
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Tho are numerous and are quite difficult to calculate. How- 
ever since the Green’s function was derived by considering only the 
fundamental equations and the boundary condition and hence is complete- 
ly independent of the proscribed blade loading or the blade ahapo, 
m 

these values 1^ can bo used for any axial flow iroblcm for which 

m 

the hub- tip ratio is tho ecoe ae that for which the 1^ are calcu- 
lated. 

Considerable simplification of the problem of oaloulating the 
m 

ie possible beoause of the ayncaetry of the Green's function. This 

function is symmetrical about /3 - z* a foot which reduces the number of 
m 

1^1 by almost one-half. Furthermore it posse a so a another sort of sym- 
metry in /3 and z in that the influence of the vortiolty at p on tho 
velocity at z is the same aa the influence of the vortioity et z on the 
velocity at . Ibis condition again reduces the number of by al- 
most one-half. This latter "reciprocal* relation oan be used with the 
radial distances o< and r for further simplification, but the relation 
is slightly different booauae of the factor in the Green's function. 

The 1^ will of course depend on the number and sizo of tho rec- 
tangles into whioh the region is divided. The reotangles should be 
small in the vicinity of the region where the blede forces act and 
can be larger away from the bledos. According to the linearized results 
given by Marble* it appears that in most cases it is not necessary to 
carry the integration more than about five blade span lengths ) 

upstream and downstream of the blades, but this range may not be suf- 
ficient for the higher approximations and it is believed advisable to 
calculate the 1^ for a somewhat greater range. 
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Ihe expression for u* is beat expressed in matrix form for mech- 
anical oaloulatlons. vre define 
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Using this notation we can write 



r 

/ 



mi 



i * ', 2 , 3 J 4- 






The four column matrices ^O^'jnuat be evaluated for eaoh iter- 
ation. The four sq.uaro matrices can bo determined once and for 
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all as soon as the subdivision of the field into aziall rootanclos 
is decided upon. 
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d) Evaluation of the IntoKral for u** 

The integi*al u ' 1 is the regular part of the integration over 
the small region containing the point at which the velocity io to be 
calculated 



Zy r + r 
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V«'e have from Eq. 33 the value of the integrand at the center 
point of the region, i .e. , at «< ( /? = r»z. as 



+ r ,n 



/ - c os Z TT 



r-n 
r* - r,_ 
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It oan be shown that this point is a saddle point of the logarithmio 
function in the integrand since it represents a maximum with respect 
to o< and a minimum with respect to (i • 

It follows then that if the region is sufficiently small and if 
the function f( °< % p>) is smooth we can use the value of the logar- 
ithmic term at tho contor as an approximate mean value (constant) so 

that an approximate value for u* • is 

z-r e- r + S' 

u" = S- In J f p) d-< dp (48) 

z-fc r-<T 

Since this integral u" represents the effect of only one small 
region the above expression is sufficiently acourate unless the 
change of vortioity in this region is much greater than anywhere else 



in the flow field 
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£© can write the integral et>ove for the /nth rectangle as 



= J m D 



( 49 ) 



where 
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Is evaluated at the /nth point and 

z-t 6 r+f 

O” - f f /?-, ffj d* d(3 (5D 

Z - 1 7- -f 

is taken over the /nth reotangle. 

The J 111 are determined once and for all as soon as the subdivisions 
of the field are deoided upon. The dP must be determined for eaoh 
iteration. 



) E valuation of the Integral for u*** 

The integral u”' is the singular part of the integration over 
the small region containing the point at which the velocity is being 
calculated* 
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«e define a new function momentarily as 
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and write the integral as 

z+e r / J' 

u"‘ — f f 9 (-<, 0 . r ) 

2 -(■ S-S 

The function g( <=*,/* »r) is regular and possesses higher order 
derivative so that it can be expanded about the oenter of the reotangle 
in a ilaclaurin's series in the two variables °< * f 3 . 
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If we substitute this expansion of g(<=*#/3) in the integral u 1 ' 
and negleot terms of fourth and higher orders we have 



= 9(0 d K, - t(0-)^ Kz + H&lt K* 



( 55 ) 
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where 
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'Ihe integrals involving the odd powers of (<*. -r ) or ( Z 3 -z) 
vanish beoause of the symmetry of the logarithmic term. The above 
integrals can be ovaluatod by straightforward integration by parts 
to givo 
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The throe constants, K^, .v,, Kj, are determined once and for all 
as soon as the subdivisions of the field are decided upon. 

By straightforward differentiation we can express the deriva- 
tives of g( o< , j ,r) in terms of % p ) and we can write the integral 
u* " for the velocity at the rectangle m as 
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d) The Complete Integral for u 

Ke have evaluated the integral for the radial velocity in thro© 
parts as given by Sqs. 42 , 49 , and 59 • Combining these results we 
can write a complete expression for the radial velocity at the center 
r,z of the rectangular subdivision m. 

= Z- [c n , I „7 +■ Cnz C z + C„,C+C n +C 

nfnn 

+ D™J + 2. 

* = C 2 ; 3 

The constants C, D and 2 are evaluated from the unhomogeneous 
part f(°<,/3) of the differential equation for the radial velooity 
(Sqs. 11 and 27) and ore to be determined in the manner described 
in the preceding soctions. These constants must bo determined for 
each iteration. 

The constants I, J and X depend only on the Green's function 
(i.e., the boundary conditions), the hub-tip ratio, and the size of the 
subdivisions and can therefore be evaluated once and for all. They are 
used in the same manner in eaoh iteration step. 
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v. Discussion 

A method of solving the hydrodynamic equations for the incom- 
pressible flov of an invlscid fluid through an axial flow turbomachine 
has been developed as an iteration procedure. The next logical step 
of this approach to the problem is the evaluation of the invariant 
terms of the iteration expression for the radial velocity. Once this 
is done for the several likely hub- tip ratios it will be relatively 
easy to calculate, by mechanical means, solutions for any axial flow 
machine with any proscribed blade. loading. 

In determining the degree of accuracy that will be required the 
restrictions (incompressibility, etc.) which were imposed to simplify 
the mathematical problem must be considered end ovaluated. The solu- 
tion should be consistent with these restrictions and of sufficient 
accuracy to indicate the proper trends of the variables. It is not 
certain that the linearized solution nests those requirenents in all 
cases. The iteration procedure con be used to determine the accuracy 
of the linearized solution and if necessary to obtain solutions of 
greater accuracy. 
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VI APP2JDIX 

A Brief Con^rlaon of the Axial Flov and Mixed Flow Problems 

The extension of the nore recent axial flow solutions to apply 
to nixed, flow machines would be extrenely difficult and would require 
consideration of several points, not significant in axial flow, which 
are of utmost importance in nixed flow. Three essential differences 
between the two problems are pointed out here in order to indicate 
which of the assumptions used in the axial flow analysis would not 
be applicable to a nixed flow analysis. 

(1) Tho nixed flow machine contains continuous wanes as con- 
trasted with tho rotor and stator blade rows in an axial 
flov machine. The vanes therefore cannot be "twisted" 
without introducing excessive tilt away fron radial so 
that the centrifugal forces cause large bending nooents 
in tho vane and prevent operation of the machine at ex- 
trenely high speeds. A geonetrlc relation oust hold be- 
tween the relative velocity, the vane forces on the 
fluid and the shape of the vane3 throughout the region 
where the vanes are present, whereas for blade rows as 
in axial flow this relation holds only in the region of 
the narrow blade and, in fact, my be concentrated in a 
"lifting line" for a good apnroxination. 



( 2 ) The larger radial velocities which naturally occur in 
nixed flow prohibit the simplifying assumption used 
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for axial flow th.it the r?dial velocity la very small 
corap, red to the axi^l velocity. 

(3) The boundary conditions of the mathematical problem 

of the nixed flow compressor are nuch nore complicated 
than those for axial flow end greatly increase the dif- 
ficulty of obtaining a solution. The variables are not 
separable in thie case. It should be noted too that 
the boundaries will probably be very different for each 
mixed flow problem whereas they are always essentially 
the sane for axial flow. 

Several type# of vnnes are possible but two special vane shapes 
are likely to be of interest. 

The first of these might be called "radial vanes". These vane# 
are generated by radial lines through the axis of rotation. Hera the 
radial force is zero. Radial vanes are necessary in a high speed 
machine because of the high centrifugsl forces. If radial velocities 
and (or) pitch angles are largo as in nixed flow then the angles be- 
tween the rane9 rnd the hub or shroud will be acute, thus increasing 
boundary layer effects. 

The second typo might be called "normal vanes". These vanes are 
generated by linos through the axis but tilted in a meridional plane 
so as to be normal to the meridional trace of the streamline at all 
points. Here the vane force has two components, one acting along the 
meridional trace of the streansurface end one aoting tangentially. 

This vane is not structurally adequate for extremely high speed rotors 
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but does hpTo the good feature th« t angles between the vane and the 
hub or shroud ro right angles, a f'->ct vr'.ich nay minimize boundary 
layer effects. It is interesting to note that one force component 
accelerates tho fluid in its path in a meridional plane and the other 
component accelerates it tangentially so that it appears that no 
"vasted* forces are present. For the axial flow machine radial vanea 
and normal vanes are about tho sane. 

The selection of either radiul or normal vanes will lead to 
great simplification of tho mathematical problem. 



V 
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b) Sona Jemecuonces of the Hydrodynamic Equations as An-plied Jo 
T urbo ray, chi no s 

The Suler equations of notion are written in vector fora and in 
c/lindrical component fora for the isentropic flow of sn inviscid 
ruid incompressible fluid acted upon by non- conservative body forces. 
The flow i3 symmetrical about the axl3. 
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By use of a vector identity the-e equations nay be written in 
^ terms of the vorticity 



1/ * jn - r(-p-) - F 



(3) 



where, since the flow i3 isentropic and adiabatic, 

? ( pj = 7 

is the gradient of the total energy of the fluid. In component form 
7 -q. (3) is 

( £. I _ /r 

t-'J" — — Z> r ' /* J r r 

W f — “ _ /^€> ^ ^ 

a. - sir£ - UTi f ^ z. 
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oeverul interesting relations between velocity, vorticity, total 
energy, vane shape, etc., can bo dorivod from these equations if we 
consider the flow between tvo concentric bounding surfaces acted upon 
by vrnes or blcde3 rot ting with an angulcr volocity . Vo will 
essune an infinite number of vanes (or blades) so that tho forces of 
the vane on the fluid nay be represented by body forces symmetrical 
about tho axi3. Furthermore since the fluid is inviocid these forces 
are normal to the vane and hence to the relative velocity so that 



If we t-^ko tho scalar product of the velocity and both sides 
of Eq. 3 these results 




(5) 



o 




or 




( 6 ) 



From the second equation of notion 





( 7 ) 



From Eqs. 6 and 7 vo see that 




, where rotating vanes act 



, where no rotating vanes act 
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If ve oon aider tho case where the fluid has constant energy 
fax upstream and is free cf vorticity far upstroan and make the stip- 
ulation that the total energy of the fluid is changed only by the 
action of the moving vanes (i.e., adiabatic flow) these 2q. 8 cen be 
integrated along each ctreorngurfoce with tho result that 






far upstream 
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Tor this case of uniform energy and no vorticity far upstream, ^ 
and cDry are constant far upstream so that: 



*(?) = ^ 



; upstream of the vanes 
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~v(u>rv) • in regions of and downstream of the 

vanes if no forces such as stationary 
vane forces occur 



If any forces, such as those resulting from stationary vanes or 



stators, have tangential components then we can only say that 
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along the stroaasurface since in this case d)d-(rv)=o. b-u-t — (ft/) -fc O 

r ' 

Oombining Bqs. 3 end 9 we have 

px JC2 - ? (edrr) - F 

and if ve multiply both sides by the vorticity 

X2 • V * S2 = 12 - V (u>rr) - _Q • W = O 

we see that 



O. • jv{iX>r v) - fJ & 
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Using only the definition of vorticity it is easily 3 hovn that 



12 ■ V ( t o r v) ~ o 



(ID 



and hence 
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and further froa 7n. 10 

J? • F ~ O 

Using the definition of vorticity wl 3q. 9 it is seen that 



For thi3 special caso wo can draw the following interesting con- 
clusions: 

Upstrean o' the Vanes i 

1. The vorticity is zero. 

2. The total energy of the fluid is unifora. 

In the Region '.here Vane Forces Acti 

1. The vorticity vector is tangent to the vane surface. (Eqs. 5 
and 13). 

2. The vorticity vector is tangent to a surface of constant 
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Conbining Uq. 14 and Uq. 3 




(16) 



total energy. (Sq. 12) 
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c. 1 -0 uouxul vorticity is tangent to the 
of a vi-xie ?nd r- concentric surface of 
c./nsireu •. of the V 

1. _..c vorticity v-ctor is t ngont to th 
vector (r^l tivo to t'.e rot tin*; rcto 
. . "'e vorticity vector is tangent to a 
total energy, this surface being the 



line of intersection 
Constant tot.,1 energy. 

o relative velocity 

v* . 

- ✓ • 

surface of constant 
strearasurfr-co. 
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